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Abstract. The paper studies stochastic integration with respect to Gaussian pro- 
cesses and fields. It is more convenient to work with a field than a process: by 
definition, a field is a collection of stochastic integrals for a class of deterministic 
integrands. The problem is then to extend the definition to random integrands. An 
orthogonal decomposition of the chaos space of the random field, combined with 
the Wick product, leads to the Ito-Skorokhod integral, and provides an efiicient 
tool to study the integral, both analytically and numerically. For a Gaussian pro- 
cess, a natural definition of the integral follows from a canonical correspondence 
between random processes and a special class of random fields. Some examples of 
the corresponding stochastic difl'erential equations are also considered. 



1. Introduction 

While stochastic integral with respect to a standard Brownian motion is a well-studied 
object, integration with respect to other Gaussian processes is currently an area of 
active research, and the fractional Brownian motion is receiving most of the attention 
[D El in El El 13 Eni [ISl etc.] The objective of this paper is to define and investigate 
stochastic integrals with respect to arbitrary Gaussian processes and fields using chaos 
expansion. 

A generalized Gaussian field X over a Hilbert space H is a continuous linear mapping 
/ 1-^ X(/) from H to the space of Gaussian random variables. The correspond- 
ing chaos space is the Hilbert space of square integrable random variables that 
are measurable with respect to the sigma-algebra generated by / G H. The 

chaos expansion is an orthogonal decomposition of H^: given an orthonormal basis 
{Cm, m > 1} in M.X, a square integrable H-valued random variable rj has the chaos 
expansion r] = X]m>i Vm^m, with r]^ = ^{rjCm) G H. 

The definition of a generalized Gaussian field X already provides the stochastic inte- 
gral X(/) for non-random / G H. As a result, given the chaos expansion of a random 
element rj from H^, the definition of the stochastic integral X(?7) requires an extension 
of the linearity property of X to linear combinations with random coefficients. An 
extensions of this property using the Wick product lead to the Ito-Skorokhod integral 
X^(/); see Definition 14.11 below. Under some conditions, the integral coincides with 
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the divergence operator (the adjoint of the Malhavin derivative) on the chaos space 

Even for non-random / there are often several ways of computing X{f). It is most 
convenient to work with a white noise over H, that is, a zero-mean generahze Gaussian 
field such that E(X(/)X((7)) = {f,g)n for all f,g eH. It turns out that, for every 
zero-mean Gaussian field X over H, there exists a different (usually larger) Hilbert 
space H' such that X is a white noise over H'. Moreover, the space H' is uniquely 
determined by X. On the other hand, every zero-mean Gaussian field X over H can 
be written in the form X(/) = ?B(/C*/), / G H, where JC* is a bounded linear operator 
on H and 05 is a white noise over H, although this white noise representation of X 
is not necessarily unique. Thus, different white noise representations of X lead to 
different formulas for computing X(/), and the chaos expansion is an efficient way for 
deriving those formulas. In particular, for both deterministic and random /, chaos 
expansion provides an explicit formula for X(/) in terms of the Fourier coefficients of 
the integrand /. 

To define stochastic integral with respect to a Gaussian process X = X{t), t G 
[0,T], we construct a Hilbert space Hx and a white noise 03 over Hx such that 
X{t) = ^{xt), where Xt is the indicator function of the interval [0,t]. The space 
Hx is uniquely determined by X; for example, the Wiener process on (0, T) has 
Hx = -^2((0,r)). Then the equality 

(1.1) r f{s)odX{s) = ^%f), fem^, 

Jo 

defines the stochastic integral with respect to X. 

In some situations, given a Gaussian process X = X(t), t G [0,T], it is possible to 
find a generalized Gaussian field X over a Hilbert space H so that X(t) = X(xt). 
Even though X is not necessarily a white noise over H, the resulting definition of the 
stochastic integral, 

[ /(t)odX(t)=X^(/), 
Jo 

coincides with the f ll.ll) . while the space H can be more convenient for computations 
than the space Hx- For example, fractional Brownian motion with the Hurst pa- 
rameter bigger then 1/2 has a rather complicated space Hx, but can be represented 
using a generalized Gaussian field over H = L2((0,T)). 

The paper is organized as follows. Section [2] provides background on generalized 
Gaussian fields, the chaos expansion, and the Wick product. Section [3] establishes 
a connection between Gaussian processes and fields. Section H] investigates the Ito- 
Skorokhod stochastic integral. Section [5] studies the corresponding differential equa- 
tions that admit a closed form solution. Section E] studies more general stochastic 
evolution equations and establishes the corresponding stochastic parabolicity condi- 
tions. 
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The main contributions of the paper are: 

(1) A connection between generahzed Gaussian fields over L2((0,T)) and pro- 
cesses that are representable in the form K{t, s)dW{s) (Theorem 13. ip : 

(2) Chaos expansions of the Ito-Skorokhod integral (Theorem 14.21) : 

(3) Investigation of the equation u{t) = 1 + u{s)odX{s) for a class of Gaussian 
random processes X (Theorem 15. ip . 

(4) A generalization of the stochastic parabolicity condition (Theorem 16. 3p . 

In particular, we establish the following result. 

Theorem 1.1. Let X be a zero-mean generalized Gaussian field over L2{{0,T)) and 
X{t) = X{xt). Then 

(a) The solution of the Ito-Skorokhod equation 

u{t) = 1+1 u{s)odX{s) 



is unique in the class of square integrable J-'^ -measurable processes and is given by 

u{t) = e^W-l^^'W. 

(b ) The Ito-Skorokhod partial differential equation 

du{t,x)=a / Uxxis,x)ds + aux{t,x) o dX{t) 
Jo 

is well-posed in L2{^1] L2(M)) if and only if 

2 

at > —EX^it). 
- 2 ^ ^ 



2. Generalized Gaussian Fields: A Background 

Let {Q, JF, P) be a probability space and V, a linear topological space over the real 
numbers M. Everywhere in this paper, we assume that the probability space is rich 
enough to support all the random elements we might need. 

Definition 2.1. (a) A generalized random field over\ is a mapping 
X : X V — s> M with the following properties: 

(1) For every / G V, X(/) = X(-, /) is a random variable; 

(2) For every a,P eR and f,g e V, X(a/ + Pg) = aX{f) + (3X{g); 

(3) // lim fn = fin the topology ofV, then lim X{fn) = X{f) in probability. 

n— >oo n— >oo 

(b) A generalized random field X is called 

• zero-mean, if EX(/) = for all f G V; 

• Gaussian, if the random variable X{f) is Gaussian for every / G V. 
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For Example, ii W = W(t), < t < T, is a standard Brownian motion on 
JF, P), then = f{t)dW{t) is a zero-mean generalized Gaussian field over 
L2((0,T)); note that 

(2.1) E|X(/„) - X(/)p = r - f{t)\'dt. 

Jo 

More generally, if is a bounded linear operator on L2((0,T)), then 

(2.2) X(/)= f {Mf){t)dW{t) 

Jo 

is a zero-mean generalized Gaussian field over L2((0,T)). In fact, by Theorem 12.6( b) 
below, every zero-mean generalized Gaussian field over Iv2((0, T)) can be represented 
in the form (12. 2p with suitable M. and W . We will also see that the fractional 
Brownian motion on [0,T] with Hurst parameter bigger than 1/2 can be interpreted 
as a zero-mean generalized Gaussian field over L2((0,T)). 

Let H be a real Hilbert space with inner product (■, ■)h and norm || ■ ||h = \/ (■, Oh- 
The following result is a direct consequence of the Riesz Representation Theorem. 

Theorem 2.2. For every zero-mean generalized Gaussian field X over a Hilbert space 
H, there exists a unique bounded linear self-adjoint operator TZ on H such that 

(2.3) E{X{f)Xig)) = (7^/, g)n, f,geU. 

Definition 2.3. (a) The operator TZ from Theorem \2.2\ is called the covariance 
operator of X. (b) The field X is ca/Zec? non-degenerate ifTZ is one-to-one. (c) A 
white noise over H is a zero-mean generalized Gaussian field with the covariance 
operator equal to the identity operator. 

Standard arguments from functional analysis lead to the following result. 

Theorem 2.4. (a) For every zero-mean generalized Gaussian field X over a Hilbert 
space H, there exist a bounded linear operator /C on H and a white noise ?B over H 
so that /C/C* is the covariance operator of X and, for every / G H, 

(2.4) Xif) = 5B(/C*/); 
as usual, K* denotes the adjoint of fC. 

(b) For every zero-mean non- degenerate generalized Gaussian field X over a Hilbert 
space H, there exists a Hilbert space Ht^ such that H is continuously embedded into 
tl-ji and X extends to a white noise over Ht^. 

Remark 2.5. (a) If X is non- degenerate and TZ : tl ^ tl is onto, then TZ has 
a bounded inverse and tl-ji = H. (b) If kevTZ is non-trivial, then we can define 
Hji as the closure of the factor space H/ker(7?-) with respect to the inner product 
{f,9)nu = (J^f,9)H, where f is the equivalence class of f in H/ker(7?.). Direct 
computations show that the generalized random field 03 over H/ker(7^), defined by 

^(7) = X(/), / G H, 

extends to a white noise over H-r,. 
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We will now discuss several connections between generalized Gaussian fields and 
Gaussian processes. 

Denote by Xt = Xt{s) the indicator function of the interval [0,t]: 

1, 0<s<t: 



(2.5) xtis) 



0, otherwise. 



With this definition, Xt2{^) ~ Xti{s) is the indicator function of the interval (^1,^2], 
t2 > h. 

Theorem 2.6. (a) If^ is a white noise over L2{M.) , then B{t) = ^{xt) is a standard 
Brownian motion on {Q, J-', P) and, for every f G ^2(1^), we have 



(2.6) 53(/) = / f{s)dB{s). 



(b) For every zero-mean non-degenerate generalized Gaussian field X over L2{{0, T)), 
there exist a bounded linear operator /C* on L2(M) and a standard Brownian motion 
W = W{t) such that, for every f G L2{M.), 

(2.7) X(/) = / {IC*f){s)dW{s). 

Jr 

Proof, (a) Direct computations. 

(b) This follows from part (a) and from Theorem 12.41 □ 

Remark 2.7. While beyond the scope of this paper, a similar result is true for multi- 
parameter processes as well. For example, if X is a generalized Gaussian field over 
L2(M^), then the same arguments show that (12. 7p holds with a Brownian sheet W . 

From now on, we assume that the space H is separable, the field X is 
non-degenerate, and JF = JF-^, the sigma-algebra generated by the random 
variables X(/), / G H. 

Definition 2.8. (a) The chaos space generated by X is the collection of all square- 
integrable random variables on [fl, JF, P). This chaos space will be denoted by M.x- 

(b) The first chaos space generated by X is the sub-space ofMx, consisting of the 
random variables X{f), / G H. The first chaos space will be denoted by M.^^\ 

It follows that M.X is a Hilbert space with inner product (^,'/7)h^ = ^{^v)y is 
a Hilbert sub-space of Mx- Moreover, the space is separable: if {/i, /2, . . .} is a 
dense countable set in H, then the collection of all finite linear combinations of X{fi) 
with rational coefficients is a dense countable set in H^'*. 

Our next objective is to show how an orthonormal basis in H^'' leads to an orthonor- 
mal basis in H^. We will need some additional constructions. 

For an integer n > 0, the n-th Hermite polynomial if„ = Hn{t) is defined by 

(2.8) H^{t) = (-l)VV2l-e-*V2. 
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Next, denote by X the collection of multi-indices, that is, sequences a = {ak, k > 
1} = {ai, a2, . . .} with the following properties: 

• each ak is a non-negative integer: G {0, 1,2,...}. 

oo 

• only finitely many of are non-zero: \a\ := ^ < oo. 

k=l 

The set X is countable, being a countable union of countable sets. By e„ we denote the 
multi-index a = {ak, ^ > 1} with ak = I ii n = k and ak = otherwise. For a & J', 
we will use the notation a\ := ai! ■ ■ ■ • Let {^1,(^25 • • •} be an ordered countable 
collection of random variables. For a G X define random variables C,a as follows: 

(2.9) ea=n^"'^^'^ 



k>l v«fc! 

where H^f. is a^-th Hermite polynomial (12.81) . 

Theorem 2.9. Let {^1,^2, • • •} be an orthonormal basis in H^''. Then the collection 
S = {^o, a G X} zs an orthonormal basis in Hx-' for every rj G Mx we have 



Proof. See [13i Theorem 2.1]. □ 

Corollary 2.10. Let ^ be a white noise over a separable Hilbert space H and let 
{mi, m2, . . .} be an orthonormal basis in H. Then {^^ = Q3(mfc), k > 1} is an 
orthonormal basis in H^^'' and, for every / G H, 

00 

(2.10) «(/) = J](/,mfc)H^(mfc). 

fc=i 



Proof. Note that E(,^fc,^„) = E(^(mfe)^(m„)) = (mfc,m„)H, so the system {S^k, k > 
1} is orthonormal in H^'' if and only if {nik, A; > 1} is orthonormal in H. If ^ G H^'', 
then ^ = ?B(/) for some / G H. By assumption, / = ^k)-R "^fc, which implies 

(I2.10p and completes the proof. □ 



If H = L2((0,T)) and / = Xt, then (12.101) becomes a familiar representation of the 
standard Brownian motion on [0,r]: 

(2.11) W{t) = £ mk{s)d^ (^j^ mk{s)dW{s)^ . 

Now, let X be a zero-mean generalized Gaussian field over a separable Hilbert space 
H. By (I2.10p and Theorem 12.4( a). we can take a white noise representation of X, 
X(/) = *B(/C*/), and get an expansion of X(/) using an orthonormal basis in H: 

00 

(2.12) X(/) = Y,^lC*f,mk)ii'B{mk). 

k=l 
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Alternatively, by Theorem 12.4( b) X is a white noise over the space Ht^ corresponding 
to the covariance operator TZ of X. If {mfc, /c > 1} is an orthonormal basis in H-/^, 
then we have an equivalent expansion of X(/): 



(2.13) :£(/) = 5^ (7^/,m,)HX(mfc) 



k=l 



We conclude the section with a brief discussion of the Wick product, as we will need 
this product to define X(/) for random /. For more details, see [8]. 

The Wick product of two arbitrary elements of Hx can be computed by the formula 



where a + P = {a^ + Pk, k > 1} and a\ = Y[k>i '^fc' ~ ai!a2'«3' ■ ■ ■ • In general, there 
is no guarantee that, for ^, G Hx, the Wick product ^ o 77 belongs to Mx- 

Similar to ordinary powers, we define Wick powers of a random variable rj G Mx'- 
^on _ fj o ■ ■ ■ o rj. Replacing ordinary powers with Wick powers in a Taylor series 
for a function / leads to the notion of a Wick function For example, the Wick 
exponential e^^ is defined by 



(2.15) e^ = J2 



n=l 



and satisfies e^'-^"'"''^ = e*^ oe*''. If 77 G H^'', then direct computations show that 

(2.16) gor, _ gr,-iEr,2_ 

3. Connection Between Processes and Fields 

Given a zero-mean generalized Gaussian field X over L2((0,T)), we define its 
associated process X{t), t > 0, by 

(3.1) Xit) = X{xt). 

Clearly, X{t) is a Gaussian process. Let /C* be the operator from Theorem 12.61 and 
define the kernel function Kx = Kx{t, s) by 

(3.2) Kxit,s) = ilC*Xt)is). 
It then follows from (12. 7p that 

(3.3) X{t) = [ Kx{t,s)dW{s) 

Jo 

for some standard Brownian motion W. Let us emphasize that, while every kernel 
K{t, s) with minimal integrability properties can define a Gaussian process according 
to (13.31) . only a process associated with a generalized field over L2((0, T)) has a kernel 
defined according to (13. 2p . where /C* is a bounded operator on L2((0,T)). Recall 
that the definition of a generalized field (Definition 12.11) includes a certain continuity 
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property, and this property translates into addition structure of the kernel function 
in the representation of the associated process. 

Now assume that we are given a Gaussian process X{t) defined by (13. 3p with some 
kernel Kx{t, s). We are not assuming that Kx has the form (13. 2p . In what follows, we 
discuss sufficient conditions on Kx{t, s) ensuring that X{t) is the associated process 
of a generalized Gaussian field X over L2((0,T)), that is, representation (13. 2p does 
indeed hold with some bounded linear operator /C* on L2((0,T)). For that, we need 
to recover the operator JC* from the kernel Kx{t, s). By linearity, if (13. 2p holds and if 
So < si < . . . < Sat are points in [0, T] and 

N-l 

(3.4) fis) = Y.'''^(Xs,^^(')-XsM) 

k=0 

is a step function, then 

N-l 

(3.5) K,*f{s) = ctfc {Kx{sk+i, s) - Kx{sk, s)) . 

To extend (13.50 to continuous functions /, the kernel Kx(t, s) must have bounded 
variation as a function of t; if this is indeed the case, then (13.50 implies that, for every 
smooth compactly supported function / on [0,T], 

(3.6) /C7(s)= f f{t)Kx{dt,s). 

Jo 

It now follows that if the partial derivative dKx{t, s) / dt exists and is square integrable 
over [0, T] x [0, T], then /C*, as defined by (13.61) . extends to a bounded linear operator 
on L2((0,T)). 

Let us now assume that the process X(t) define by (13. 3p is non-anticipating, i.e. 
adapted to the filtration {J-"^ , < t < T} generated by the Brownian motion W{s). 
Then Kx{t, s) = for s > t and (13. 3p becomes 

(3.7) X{t)= [ Kx{t,s)dW{s). 

Jo 

Note that in this case we have 

/>min{t,s) 

(3.8) E(X(t)X(s))= / Kxit,T)Kxis,T)dT. 

Jo 

For such processes, formula (13.50 and the conditions for the continuity of the corre- 
sponding operator /C* must be modified as follows. 

Theorem 3.1. Assume that the process X{t) defined by liS. 3\) is n on- anticipating. 



(a) If f is a step function (3^, then 



N-l 



(3.9) 



^*/(s) = Yl U^>+i(^) - XsA^)) (aiKx{si+i, s) 

1=0 
N-l 

+ ^ ak{Kxisk+i,s) - Kxisk,s) 



k=i+l 
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(b) If the function Kx{-, s) has bounded variation for every s and lim Kx{s+6, s) = 

(5— >0, S>0 

Kx{s~^, s) exists for all s G (0,T), then 

(3.10) ]C*f{s) = Kx{s+, s)f{s) + fmxidt, s) 

for every continuous on [0,T] function f. 

(c) If the function Kx{t, s) has the following properties 

(1) Kx is continuous and non-negative for < s < t < T , and sup Kx{t,t) < 

0<t<T 

(2) For every fixed Sq, the function Kx{t, Sq) is monotone as a function of t and 
the partial derivative K^^\t, s) = dKxit, s)/dt exists for all < s < t < T ; 

(3) There exists a number Ki = Ki{T) such that 

(3.11) sup / Kx{T,s)\K^^\t,s)\ds < Kf, 

0<t<T Jo 

then the corresponding operator IC* defined by equation (13.101) is bounded on L2{{0, T)) 
and the operator norm ||/C*|| o//C* satisfies 

(3.12) \\IC*\\^ <{Ko + K,f. 

Proof, (a) By assumption, Kxit, s) = for s > t. Fix an s such that s G {sj, Sj+i] 
for some j = 0, . . . , iV — 1. By (I3.5p we have for this value of s 

N-l 

^*fis) = X] "-kiKxi-Sk+i, s) - KxiSk, s)) 

k=0 

= ajKxisj+i,s) + ^ ak{Kxisk+i,s) - Kxisk,s)). 
k=j+i 

Since Xs^+ii^) ^Xs^i^) is the indicator function of the interval (s^, Sk+i], 03. 9p follows. 

(b) Under the additional assumptions on the kernel Kx, (13.101) follows from (13.91) after 
passing to the limit max |so+i — sA 0. 

j=0,...,N-l 

(c) With no loss of generality, we can assume that Kx{s+,s) = and K^^^ > 0; 
otherwise, we replace Kxit, s) with either Kx{s~^, s)—Kx{t, s) or Kx(t, s)—Kx{s^, s). 
Let 5^ be a smooth compactly supported function on (0, T). It follows from (I3.10p that 
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Then we use the Cauchy-Schwartz inequahty and the properties of K^^^ : 







ds 



ds 



<[ [ K^^\r,s)dT I K^^\T,s)g\T)dT ds 

Jo Js Js 

< {Kx{T,s)-K^{s,s)) J K^^\T,s)g\T)dT ds 

< j^^ Q\^{T,s)K^'\t,s) ds^ g\r) dr < Kl{T)\\g\ 



2 

i2((0,T))- 



□ 



Here are several examples of processes covered by part (c) of Theorem I3.1[ 

Example 3.2. Assume that Kx{s^,s) = and K'^^\t,s) > 0. Let R{t,s) 
E{X{t)X{s)). By ([31]), 

^^^ = I^K^{T,s)K('\t,s)ds 



and therefore 

(3.13) K'f{T) = sup 



2,^^ 9RiT,t) 



0<t<T 



dt 



In particular, for the fractional Brownian motion on [0, T] with the Hurst pa- 
rameter H > 1/2, 

i?(T,t) = l(T2^ + t2^^-(r-tf^), 

and (see Nualart [HI Section 5.1.3]) has representation (13.71) with 
K^{t, s) = cJh-1] s\-" [\t - sf-h^-^ dr, 



where 

C -( ^H^il'H) 
^ \r {H + l)r{2-2H) 

and r is the Gamma-function. In this case, Kx{s+,s) = and K^^\t,s) > 0. By 

dan, 

(3. 14) (T) = 2H 2^~^^T^^-\ 

The bound Ki{T) is asymptotically optimal: if if = 1/2, which corresponds to the 
standard Brownian motion, the right-hand side of (13.141) is equal to 1. 

Example 3.3. Assume that 

K^it,s) = p{it~sr)Xtis), a>0, 

where p is a non-negative, monotone, continuously differentiable function on [0,T]. 
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If p is non-increasing, then Kq = p{0) and Kf{T) = p(0) (p(0) -p(T")). In particular, 
consider the stable Ornstein-Uhlenbeck process 

(3.15) dX{t) = -bX{t)dt + dW{t), X{0) = 0, 6 > 0, 

so that X{t) = J^e-^'^'"'>dW{s) and Kx{t,s) = e-^^'-'l For this process, ||/C*|| < 
1 + Vl — e~^^. This bound is asymptotically optimal: as 6 ^ 0, the process becomes 
W, and the upper bound on ||/C*|| becomes 1. 

If p is non-decreasing, then Kq = p(0) and Kf(T) = p{T'^)(^p{T") — p(0)). In partic- 
ular, consider the unstable Ornstein-Uhlenbeck process 

(3.16) dX{t) = bX{t)dt + dW{t), X(0) = 0, 6 > 0, 

so that X{t) = e''(*-")rfl^(s) and Kx{t,s) = e^^^-'\ For this process, ||/C*|| < 

1 + lye^^X^^^^"^^^^- This bound is asymptotically optimal: as 6 — 0, the process 
becomes W , and the upper bound on ||/C*|| becomes 1. 

The following theorem establishes a connection between a zero-mean Gaussian process 
and white noise. 

Theorem 3.4. For every zero-mean Gaussian process X = X{t), t G [0,T] with 
X(0) = and the covariance function R(t,s) = E(X(i(:)X(s)), there exist 

(1) a Hilbert space Hr containing the indicator functions Xu 

(2) a white noise 03 over 

such that X{t) = ^{xt)- 

Proof. Let the Hilbert space H/j be the closure of the set of the step functions with 
respect to the inner product 

{Xti,Xt2)HR = Rihih)- 
Define a generalized Gaussian field 03 over by setting 

(3.17) ^{xt)=X{t), 

and then extending by linearity and continuity to all of Hr. With this definition, ^ 
is a white noise over Hpi- ^ 

By analogy with (12.61) . if X is a generalized Gaussian field over a Hilbert space H 
containing Xt, t G [0, T] and X{t) is the associated process of X, then f{s)dX{s) 
can be an alternative notation for X(/). 

If X(t) is the associated process of a zero-mean non-degenerate generalized Gaussian 
field X over H = L2((0,T)), and 71 is the covariance operator of X, then R(t,s) = 
{Tixt,Xs)L2{(Q,T)) and the space coincides with H?^ from Theorem 12.41 

If R{t,s) = min(t, s), then {xti,Xt2)nR = (Xii, Xt2)L2((o,T))- That is, for the Wiener 
process, Hr = L2((0,T)). For the fractional Brownian motion, the space can 
be characterized using fractional derivative operators [13]. For a general Gaussian 
process X with covariance function R, an explicit characterization of the space Hr 
is impossible. 
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4. Stochastic Integration 

In the definition of a generalized random field X over a Hilbert space H, we consider 
random variables for non-random / G H. In this section, we define ^{f]) for 

H-valued random elements 77. 

One possible way to proceed is to write the chaos expansion of r], r] = Va^a, Va ^ 
H, and then define ^{f]) as a corresponding linear combination of X{T]a). In the case 
H = L2((0, T)), another possibility is to take a partition = to<^i<---<^iv = ^ 
of the interval [0,T] and approximate rj{t) with a sum J2iLiV{'ti){Xti — where 
t* e [tj_i,tj], and then approximate X(77) with the corresponding linear combination 
of ^iXu) -X{xu_,). 

Either way, we need to address the following question. By definition, if a, (3 are real 
numbers and /, g are elements of H, then + Pg) = + PX{g). But what if 

a and f3 are random variables? One possibility would be to keep the same linearity. 
In the case H = L2((0,T)) this would imply 

(N \ ^ 

1=1 / i=l 

where X{t) = is the associated process of X. While natural, this extension of 

the linearity property can lead to ambiguities in the definition of the corresponding 
stochastic integral. Indeed, let *B be a white noise over L2((0,T)), and let ?7(t) = 
!B(xt). By Theorem 12.61 is a standard Brownian motion W and, as we know, the 
limit of the sum W{t*){W{ti+i — W{ti)) depends on the location of the points 
f* 

Let us now consider an alternative to (14.11) : 

(TV \ ^ 

i=l J 1=1 

This time, if we take 77 (t) = X{t) = W{t), a standard Brownian motion, then the 
limit is the Ito integral W{t)dW{ty, it is equal to {W\T) - T)/2 = (I^(T))*V2 
and does not depend on the location of the points t*. 

Accordingly, we adopt the following convention: if X is a generalized Gaussian field 
over H, then, for every /, (7 G H and all random variables ^, t], 

X(e/ + w) = ^oX(/) + r/oX(^7); 

recall that, by assumption, the underlying probability space (f2,jF, P) is such that 
J-" is generated by X. With this convention, we proceed with the definition of the 
stochastic integral X(r7) for H-valued random elements r] using the chaos expansion. 
There are at least three advantages of the chaos approach over time partitioning: 

(1) generality: spaces other than L2((0,T)) can be considered; 

(2) possibility to use weighted chaos spaces, which eliminates many questions 
about convergence; 
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(3) computational efficiency: the Wick product must only be computed for the 
basis elements ^a- 

Because we are not restricted with the choice of H, we will assume that X = ?B, 
a white noise over H. Let {m^, A; > 1} be an orthonormal basis in H. Define 
= ^("^fc) and ^Q,, a G X, according to fl2.9p . 

By Theorem I2.9[ every H- valued random element rj with IE||?7||h < has chaos 
expansion 

(4.3) V = Y1 = ^(^^") ^ 

Definition 4.1. Let t] be an H-valued random element with chaos expansion (^7^. 
The Ito-Skorokhod stochastic integral of rj with respect to 03 is 

(4.4) ©^(r/) = 5^^(r/„)oe„, 
where o is the Wick product. 

Since every generalized Gaussian field and every Gaussian process can be represented 
using a white noise over a suitable Hilbert space, formula fl4.4p defines stochastic 
integral with respect to any Gaussian process or field. We will see below that this 
formula also provides a chaos expansion of the integral in terms of the chaos expansion 
of the integrand; note that (14.41) is not a chaos expansion in the sense of (14.31) . The 
two immediate question that are raised by the above definition and will be discussed 
below are (a) the convergence of the series, and (b) the dependence of the integral on 
the choice of the basis in H. 

We start by deriving the chaos expansion of the integrals without investigating the 
question of convergence. 

Theorem 4.2. Let rj be an H-valued random element with chaos expansion lli4-3\ ), 
and assume that 

oo 

(4.5) r]a = y^?7a,fcmfc. 

k=l 

Then 

(4.6) = E ( E V^^Va-e„k] U 

a£l \k=l / 

Proof. By (14.51) and linearity, keeping in mind that both rja and are non-random, 

oo oo 
k=l k=l 

Therefore, 

oo oo 

(4.7) ® (^) = E E ''"'^^^^ ^ = E E v^"*: + 1 ^".fc ! 

oGX k=l oSX k=l 
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recall that is the multi-index with the only non-zero entry, equal to one, at position 
k. By shifting the summation index, we get (14 .6 1) . Note that, for every a G X, the 
inner sum in (14.61) contains finitely many non-zero terms. 

□ 



Now, let us address the questions of convergence and independence of basis. The 
Cauchy-Schwartz inequality implies that if 

(4.8) ll^alln < oo> 

then 55* (77) e Hsg. Further examination of (14. 6 p shows that, for every rj satisfying 
]8]), 05^(77) = S{ri), where 5 is the divergence operator (adjoint of the Malliavin 



derivative), and therefore 03* (77) does not depend on any arbitrary choices, such 
as the basis in H; for details, see Nualart [13] or Watanabe [IS]. In particular, if 
H = L2((0,T)), then ^^{rj) is the Ito-Skorokhod integral of rj in the sense of the 
Malliavin calculus. On the other hand, (14.61) . if considered as a formal series, allows 
the extension of to weighted chaos spaces, similar to those considered in [111 [12] ; 
we leave this extension to an interested reader. 

Remark 4.3. It is also possible to define 

(4.9) ^°(r/) = ^^(r/J.e„, 

where ■ is the usual product. To understand the structure of this integral, note that 
the Malliavin derivative © of satisfies 



a—eic^kj 



k=l 



this follows directly from the definition of D [14, Definition 1.2.1] and the relation 
H^[x) = nHn-i{x). On the other hand, direct computations show that 

00 

ael k=l 

and, using the following property of the Hermite polynomials, 

Hi{x)Hn{x) = Hn+i{x) +nHn-iix), 

we get 



(4.10) = ( n ""^^^^^^^^^^^ = v/^e... + v^ki 



or 



(4.11) ^°(r7) = J] 5^ [V^k^la-e^k + V^^r]a+e„k) i 



a&I \k=l 
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^45 a result, we use to re-write ^.11^ as 

(4.12) Q3°(r7) = ^^ir]) + Y,iVaMa)n. 

That is, 55° is a sum of the Ito-Skorokhod integral plus the trace of the Malliavin 
derivative — a representation characteristic of the Stratonovich-type integrals [TJj. In 
particular, if ^ is a white noise over L2{{0,T)) and W{t) = ?B(xt), then ^°{W) = 
W\T) = /q W{t) odW{t). More generally, if r] = 7]{t) is in the domain of the 
Malliavin derivative, then 

(4.13) ?B°(r/) = ^%ri) + I Btrjdt, 

Jo 

where 

(oo 
^/a^ia-e^mk{t) 
k=l 

unlike the Ito-Skorokhod integral, though, condition is not enough to ensure the 
existence of^°{r]) as an element of M.x- When H is the Hilbert space of functions 
on an interval [0,T], square integrable with respect to a (not necessarily Lebesque) 
measure fi, the sufficient conditions for the Stratonovich integrability are discussed in 
[H Chpater 3]. 

Let X be a zero-mean non-degenerate generalized Gaussian field over a separable 
Hilbert space H. As we mentioned earlier, by Theorem 12.4( b) on page HI X is a white 
noise over a bigger Hilbert space H-n, and then X*(?7) can be defined using fl4.4p . If 
the space H-n is difficult to describe, one can use representation (12. 4p from Theorem 
12.4( a) and derive an equivalent formula for the stochastic integral: 

(4.14) X^(r/) = ?B^(/C*r7). 
for every (*B, H)-admissible r]. 

Unlike (14.41) . representation (14.141) is not intrinsic: the operator /C* and the white 
noise *B are not uniquely determined by X. On the other hand, in many examples, 
such as fractional Brownian motion with the Hurst parameter bigger than 1/2, it 
is possible to take H = L2((0,T)), and then (14.141) becomes more convenient than 
(14.41) . To derive the chaos expansion of X^{ri) using (I4.14p . fix an orthonormal basis 
{mk, A; > 1} in H, define C,k = ^(^fc), and consider the corresponding orthonormal 
basis {^a, a G X} in Hsg constructed according to (12. 9p . It follows from (14.61) that 

(4.15) = E ( Ev^^a-e.,fc ) ^a, 

aeX \fc=l / 

where 



If H = L2((0,T)), then (Hl5ll becomes 
(4.16) ^%v) = J2ij2^^( f 



r]a-e^,{t){ICmk){t)dtj I 
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where rjait) = K(^ri(t) ^a)- In this case, by analogy with the Brownian motion, 

Jq ri{s) o dX{s) can be an alternative notation for 3C^{ri), where X{t) is the asso- 
ciated process of X. 

5. Stochastic Evolution Equations with Closed-Form Solutions 

In this section we consider stochastic differential equations driven by a white noise ^ 
or some other zero-mean generalized Gaussian random field X over a Hilbert space H. 
To introduce time evolution into the stochastic integral, we use the function Xt, the 
indicator function of the interval [0, t], and define time-dependent stochastic integrals 

(5.1) ^tiv) ■■= ^ti^) ■■= ^'(Vxt). 

These definitions put an obvious restriction on the Hilbert space H, which we call 
Property I: H contains Xt, t G [0, T], and, for every r/ G H and every fixed t, the 
(point-wise) product rjxt is defined and belongs to H. There is a more significant 
restriction on H, which we illustrate on the following equation: 

(5.2) u{t) = 1 + ^*(n), < t < T, 

where ^ is white noise over a Hilbert space H with Property I. Let us assume that 
the solution belongs to M<s so that u{t) = J2a€:j'^ai^)^o^ each Ua is an element of 
H. By (15. ip . we can re- write (15. 2p as 

(5.3) u{t) = l + ^\uxt), 
and then (14.60 implies 

oo 

(5.4) u^{t) = l + Y, 

k=l 

Thus, the expression {ua-ekXtyn^^k)!!, as a function of t, must be an element of H, 
and the Hilbert space H must have another special property, which we call Property 
II: for every /, (7 G H, the inner product {fxt,g)ii, as a function oft, is an element 
o/H. The space Hn corresponding to a zero-mean Gaussian process with covariance 
R contains step functions by definition, but the point-wise multiplication is a more 
delicate issue, especially if is smaller than L2((0,T)). On the other hand, the 
space L2((0, T), /i), with /i((0,T)) < 00, has both Property I and Property II, which 
follows from the Cauchy-Schwartz inequality. 

Theorem 5.1. If X is a zero-mean generalized Gaussian field over L2{{0,T)), then 
the solution of the equation 

(5.5) u{t) = 1 + X^tiu) 
is unique in L2((0, T); H^) and is given by 

(5.6) M(t) = e^^W, 

where is the Wick exponential function Ii2.15\) and X{t) = X{xt) is the associated 
process ofX. 
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Proof. Let = 03 (/C*/) be a white noise representation of X over L2{{0, T)). We 
start by establishing uniqueness of solution in L2{{0,T)]M.<g), which, because of the 
inclusion C Hsg, is even stronger. By linearity, the difference Y{t) of two solutions 
of (J53]) satisfies Y{t) = If Y{t) = ^ ya{t)^a, then flCTD implies 

ft 



(5.7) = ^ a/o^ / 

k>i 



where rhk = JCruk. In particular, if \a\ = 0, then ya{t) = for all t. By induction 
on ya{t) = for all a G X: if ?/q, = for all a with |a| = n, then, since 
|a — efcl = |a| — 1, equality (15. 7p implies ?/a = for all a with |q;| = n + 1. 

To establish fl5.6l) . let 

Mfc(t)= / {ICmk){s)ds. 
Jo 

By fl2llD and I^J^, 

oo 

x(t) = ^ii4(t)efc, 

k=l 

and, because of the independence of for different k, 

k>l ael 

where 

oo 

fc=l 

Similar to fl5.7l) . we conclude that if the solution u = u{t) has the chaos expansion 
u(t) = Yl 'Ua(t)C,a, then Ua(t) = 1 if |a| = and 

(5.8) Uait) = y^ jQ^ I Ua-ei,{s)rhk{s)ds 



fc>l 



if |a| > 0. Then direct computations show that 

Ua{t) = — ^=^, |a| > 1, 



satisfies (15. 8p : 



dt v«! 1 • 

fc=l * A;=l jT^fc 



> VQfc^fc(^) = V V^mfc(t)Mc.-£fc(t). 



□ 
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Corollary 5.2. The solution of 



u(t) = uo+ / a(s)M(s)(is + (o"jt? 

^0 



k=l 

where Xk are independent generalized Gaussian random fields over L2((0,T)) and 
a G Li((0,T)), cTfc G L2((0,T)) are non-random, is 

N 



u{t) = Uoexp a{s)ds^ exp^ Xk{t) 



where Xkit) = Xk{(TkXt)- 

Theorem 15.11 is a generalization of the familia result that the geometric Brownian 
motion u{t) = e^(*)-(*/2) = g^^W satisfies u{t) = 1 + u{s)dW{sy. by (ESD and 
f l2.16p . for a class of zero- mean Gaussian processes X = X{t) with covariance function 
R{t,s), and with a suitable interpretation of the stochastic integral, the solution of 
the equation u{t) = 1 + u{s) o dX{s) is 

u{t) = e-^(*)"^(o)-|(^(*'*)-2i?{t,o)+i?(o,o))_ 

Note that the proof works without assuming that X is non-degenerate or X is non- 
anticipating. On the other hand, the special interpretation of the integral is essential: 
if X is also is a semi-martingale, then the traditional Ito integral ri{t)dX{t) can also 
be defined, and the solution of the corresponding equation U{t) = 1 + J^^ u{s)dX{s) 
is the Dolean exponential 

U{t) = e-^(*)-^(o)-3Wt. 

To see that u and U can be different, consider the stable Ornstein-Uhlenbeck process 
flXT^ . for which u{t) = exp (X(t) - (1 - e-^*)/(26)) and U{t) = exp {X{t) - t/2). 
The reason for this difference is clear: while U satisfies the equation dU{t) = 
-bX{t)U{t)dt + U{t)dW{t), u{t) satisfies du{t) = -hX{t) o u{t)dt + u{t)dW{t) (it 
is known that dW = odW; see [S]). 



As an application of Theorem 15. let us find the classical, square-integrable solution 
of the stochastic partial differential equation 

(5.9) u{t, x) = uo{x) + a / u^xis, x)ds + aX'^{ux{-, x)), t > 0, a; G M, 



with a smooth compactly supported initial condition uq and constant a > 0, a & 
Let 

u{t,y) = —= / e~'^^u{t,x)dx, i = v^^. 



/27r 

be the Fourier transform of u{t,x). Then, by linearity, 

u{t, y) = Uo{y) - ay^ / u{s, y)ds + \aXt{u{-, y)), 
Jo 

and therefore 

u{t, y) = Uq exp (-ayh + ^a'^y'^R{t, t) + \yaX{t) 
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where X{t) is the associated process of X and R(t,s) = E(X(t)X(s)). With the 
notation 

2 

r{t) = at-^R{t,t), 
the classical solution of (15.91) becomes 

u{t, x) = u{t, X — aX{t)), 

where 

In particular, we get the following parabolicity or non- explosion condition for fl5.9p : 

2 

(5.10) at>yi?(t,t). 

If X{t) = W{t), then R{t,t) = t and flHJOD becomes 

(5.11) 2a>(T2. 

If X{t) = W^{t), H > 1/2, then R{t,t) = and flSTTOj) becomes 

(5.12) t^/f-i < 2a/(T2; 
this condition also appears in [7]. 

For the stable Ornstein-Uhlenbeck process (I3.15p . i?(t,t) = (1 - e"^^*)/(26), and 
condition (I5.10p becomes 

- Aht^ ' 

which is equivalent to ( 15. lip . 

For the unstable Ornstein-Uhlenbeck process (13.161) . condition (I5.10p is 

(5.13) at > ^ (e^"* - 1) . 

If (I5TTD holds, then (157[3D holds for sufficiently small t; if (l5TT|) fails, so does (l57[3l) 
for all t > 0. 

The traditional Ito version of (15. 9p with the Ornstein-Uhlenbeck process is 

du{t, x) = {auxxit, x) ± h(yX{t)ux(t))dt + aUxdW{t), 

This equation is well-posed if and only (15.110 holds, and this condition does not 
depend on h. 

We can also write the classical, square-integrable solution of 

(5.14) u{t,x) = Uq{x) + I a{s)Uxx{s-iX)ds + X1{aUx{- -iX)), t > 0, x G M, 
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with a smooth compactly supported initial condition uq and continuous functions 
a{t), cr(t). Indeed, let 

A{t) = [ a{s)ds, X,{t) = X(axt), 
R„{t,s) = E(X,(t)X,(s)), r,(t) = A{t) - ^R^{t,t). 

Then 

u(t, x) = u{t, X — X„{t)), 

where 

The parabolicity condition for (15.141) is 

r.{t) > 0, 

which in general cannot be simplified much further: when a depends on time, there 
is no easy connection between R and R^- 



6. Chaos Solution of Stochastic Evolution Equations 

Let X^, £ > 1, be a collection of independent, zero-mean, generalized Gaussian ran- 
dom fields over L2((0,T)) on a probability space (fi,JF, P). By Theorem 12.61 on page 
[5l there is a collection {We,i > 1} of independent Wiener processes and {JCg, i > 1} 
of bounded linear operators on L2((0,T)) such that 

(6.1) Xeif) = r {lClf){t)dW,{t). 







With no loss of generality, we assume that the sigma-algebra JF is 
generated by the random variables W(,{t) , £ > 1, t G [0,7]. Note that 
Xiit) = Xi{xt) is not necessarily adapted to filtration of the corresponding Wiener 
process We{t). 

Introduce the following objects: 

(1) {X,H,X'), a triple of Hilbert spaces such that X' is the dual of X relative 
to the inner product in H. To simplify the notations, we use (-,■) to 
denote both the inner product in H and the duality between X 
and X'. 

(2) A, a bounded linear operator from L2((0,r);X) to L2((0, T); X'). 

(3) Mi, i > 1, a. collection of bounded linear operators from L2((0,T);X) to 
L2((0,T);X'). 

(4) uo e L2in;H), f G L2(l^; L2((0, T); X')), ge E L^iSl; L^^i^Q^T); X')). 

(5) The Fourier cosine basis in L2((0,T)): 

(6.2) m,{s) = ^- m,{t) = ^co^[^^^^j^Yk>l- 0<t<T. 



Also define nikiit) = (/C^mfc)(t). 
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(6) J ^ a collection of multi-indices a = {aki, k,i> 1) with non-negative integer 
entries and the finite length \cx\ = '^kiCXM < oo. The two special multi- 
indices are (0) with all zero entries and with eij{ij) = 1 and e^i^ij) = 
otherwise. We also write a! = H/c^Q^fc^'- 

(7) The random variables ^ki = Jq ''nk{t)dW£{t) and 



(6.3) = n 



k,e> 



where, for an integer n > 0, Hn = Hn{t) is the n-th Hermite polynomial 
(12. 8p . By Theorem 12.91 the collection {^^t? ol G J}^ is an orthonormal basis 
in L2(fi). 

(8) the characteristic set of ex. with \cx.\ = n: Sa = {(/ci, ^i), . . . , (A;„, £„)}, 
where ki < k2 < ■ ■ ■ < kn, ii < ii+i if h = fcj+i, and 

^ ^ ^fci^i O Cfcafe O • • • O 



Consider the following stochastic evolution equation: 



(6.5) u{t) = uo 



+ [ Au{s)ds+ [ f{s)ds + J2^h{Meu + ge). 
Jo Jo 



Assume that the random element Uq and the processes u, /, have chaos expansions 
(6.6) 



Wo = Uo,c.U, fit) = ^ Uit)U, geit) = ^ gLocit)ic 

ol£J ol^J ol£J 



Substituting into (16. 5p and using fl4.16p . we conclude that the (deterministic) func- 
tions Uct{t) satisfy 

Ucx{t) = Uo^a + / {Au^{s) + fa{s))ds 

(67) ^° . 

+ 22 / i^e^a-e{k/){s) + ge,a-e{k,e){s))mke{s)ds. 

Definition 6.1. A chaos solution of equation (16.50 is a formal series u{t) = 
^^g^ MQ,(t)^ct; where the random variables o,re defined by (16. 3p and the deter- 
ministic functions Ua satisfy (16. 7p . 

Theorem 6.2. Assume that 

• for every Uq E H and F G L2((0, T); X'), the deterministic evolution equation 

(6.8) U{t) = Uo+ [ AU{s)ds + f F{s)ds 

Jo Jo 

has a unique solution U G L2((0,T);X); 
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(6.9) ess sup \fhki{t)\ < oo 

iG(0,T) 

Then equation (16.51) has a unique chaos solution and every Ua is an element of 
L2((0,r);X). 



Proof. Note that, for |a| = 0, ([62D is 

(6.10) M(o)(t) = Mo,(o) + / Au{^a){s)ds+ / f(o){s)ds; 

Jo Jo 

by assumption, equation (16.101) has a unique solution n(o) G L2((0,T);X). 

We now proceed by induction on \a\. Assume that (16. 7p has a unique solution Ua G 
L2((0,T);X) for all a with \a\ < n. Then, for a with \a\ = n + 1 we have \(x — 
e{k,i)\ = n so that MiUa~e{k,e)'f^ke G L2((0, T); X'). With only finitely many terms 
in the sum on the right-hand side of (16. 7p . the assumptions of the theorem now imply 
that (16.71) has a unique solution Mq, G L2{{0, T); X) for all a with |q;| = n + 1. □ 



In general, condition ( 16. 9p is necessary: there is no guarantee that ffhke G 
L2((0,T);X') for every / G L2((0, T); X'). More information about the operator 
A makes it possible to remove this condition; see Theorem 16.31 below. 

Assume that Xi{xt) = Jq Ki{t, s)dW£{s). Then a sufficient condition for (16.91) to hold 
is 

(6.11) ess sup lKi{t,t) + / 

te{o,T) V Jo 

this is the case for the fractional Brownian motion with H > 1/2 and for the Ornstein- 
Uhlenbeck process (stable or unstable). 

A theory of solutions of (16. 5p in weighted chaos spaces can be developed in complete 
analogy with [HI [12] ; we leave this development to an interested reader. Instead, we 
will investigate when the chaos solution is in fact square-integrable, that is, when no 
weights are necessary. This investigation will lead to an extension of condition (I5.10p 
on page [19] to equations with variable coefficients. 

Denote by the norm of the operator ICi in L2((0,T)). 
Theorem 6.3. Assume that 

(1) The initial condition Uq and the processes f,g£ are deterministic and 

Yl I ^]\\9i{t)\\%dt < ^■ 

(2) There exist positive numbers 5 a and Ca such that, for allv E X and t G [0, T], 

(6.12) \\A{t)v\\x' <Ca\\v\\x, {A{t)v,v) + 5A\\v\\\<CA\\v\\l. 



dKejt, 
dt 



ds ] < oo; 



STOCHASTIC INTEGRALS AND EVOLUTION EQUATIONS 23 

(3) There exist a non-negative number 5q < 5a and a positive number Cq such 
that, for all V E X and t G [0, T], 

(6.13) 2{A{t)v,v) + Y,^l\\Mt{t)v\\l + 54v\\\<C4v\\]j. 

i>i 

Then the chaos solution of fl6.5p satisfies 
(6.14) 

sup nnml + ^o [ E\\uit)\\j,dt<CiCA,6A,Co,T)(\\uo\\l+ [ Umx'dt 
0<t<T Jo ^ Jo 



Proof. By assumptions of the theorem, fl6.7p takes the form 

U(Q)(t) = Uo+ / Au(Q){s)ds + / f{s)ds, |q;| = 0; 
Jo Jo 

(6.15) '^e{ij){t) = Aue{.i,j){s)ds + (^MjM(o)(s) + gj{s)^mij{s)ds, \ol\ = 1; 

Ua{t) = / Auc{s)ds+ ^ ^/a^e / Mi,Uoc-e{ki){s)mki{s)ds, \ol\ > 1. 

By f l6.12p . the operator A generates a semi-group ^t,s, and the solution of (16. 8p is 

?7(t) = <l>i,of/o + f <^t,sF{s)ds. 
Jo 

By induction on \cx\ we conclude that if |q:| = n, {(fci, ii), . . . , {k^, in)} is the charac- 
teristic set of a, and Vn is the set of all permutations of {1,2,..., n}, then 

M«(t)=^5Z/ / •••/ '^MnM,^(„)---'^'s2,si(M,^(^,M(o)(si) 

(6.16) val^g^^Jo ^0 ^0 ^ 

+ 5'<?,(i)(si)j(/C£^(^jmfc^^^^j)(s„) ■ ■ ■ (/C£^(i)mfc_^(i) )(si)rfs", 
where (is" = dsi . . . dSn- We then re-write (I6.16P as 

(6.17) u^{t)= f G{t,&^;s^''^)m^{s^''^)ds'', 

J[0,T]" 

where 



+ 9h (^^(1) )) Xs,(2) (s<t(1) ) ■ ■ ■ (Sa(n) ) , 

(6.19) m„(s(")) = J— V (/C^ /ifcj(s^(i)) ■ ■ ■ (/C£„/ifcJ(s<,(„)). 

yoilnl — 



(6.18 



and 
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From f l6.17p and the definition of tlie function G, we conclude tliat 



(6.20) 



Similarly, 



(6.21) 



\cx\=n 



£i,.../„=i \i=i 




^0 



S2 



^ts Ml ■ ■ ■ $ 



S2,Sl 



M^iM(o)(si) +gi-,isi] 



H 



ft ^ I " \ ft fS pSn rS2 

E / iKwiii*< E / / / ■■■/ 

|c«|=n 



$s,.„M£„ ■ ■ (m^,M(o)(Si) +^£,(Si)) 



ds^ds. 



X 



For n > 1, denote by F^{t) and F^(t) the right-hand sides of flOOj) and flOTj) . 
respectively. For = 0, define Ff^{t) = \\u(o){t)\\'jj, F^{t) = \\u(o){s)\\xds. Then 

oo „y oo 

(6.22) nHt)\\l<J2Fn{t), / E||n(t)||^dt<E^n^m- 



n=0 



For brevity, introduce the notation 

i-T 

I 



Jo Jo 



x,dt. 



Then assumption (16.121) implies 
(6.23) 



sup F(f (t) + 5oi^o^(T) < Ci(C^, 5^, T) I. 
o<i<r 



For n > 1, we find using (16.131) that 

+ Oo — < CoF„ [t) 



dt dt 
£i,...A>i \i=i 



l-t PSn-l pS2 

Jo Jo Jo 



(6.24) 



/n+l 

E n^: 



t fSn 



fi,...A+i>i \i=i 




Jo 



l|M,„^,$i,,„M,„ . . . M,,M(o)(si) + (si) l|^c^s^ 



Then, after summation in n and integration in time, 

N N 



J2 {Fnit) + SoF^it)) <Co f J2 ^n{s)ds + C,{Ca. 6a, T) I 

n=l n=l 
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for every > 1. Applying Gronwal's inequality, we get f l6.14p from (16.221) . Theorem 
is proved. □ 



Example 6.4. Let X be the Sobolev space H^{W^) and H = L2(M^). Consider 
the following equation driven by a single fractional Brownian motion with H > 
1/2 [there should be no difficulty distinguishing between H a space and H a Hurst 
parameter] : 

du{t,x) = ^ai,(t,x)^^^^dt + ^(Ti(t,x)^^^^odW^^(t), 0<t <T. 

i,j = l i=l 

In this case, = 2H 2^-^" T^^~^: see (KW\ on page[IUl Condition fICTD becomes 

d 

(6.25) 6o\y\^ < ^ (a,,(t,x) - H2^-^''T^''-^a,{t,x)a,{t,x))y,y, < Co\y\' 

for all t G [0, T] and all x,y G M'^. Let us now compare (16.251) with (15.121) on page 



if a, a are constants and d = 1. If H = 1/2, then (16.251) becomes (15.111) . which is 
( KT2\i with H = 1/2 (recall that 6q can be zero). If i7 > 1/2, then (lOSj) becomes 
2a/cr2 > H2'^~-iHrp2H~i^ Mhich is slightly stronger than (l532|) because 1 < H2'^-'^" < 
1.07 for 1/2 < /J < 1. 

Example 6.5. Let X{t) = X{xt) and consider the equation 

(6.26) du{t, x) = aUxxdt + cru^ o dX{t), x G M, < t < T, 
with constant a, a. Condition (16.131) in this case is 

2 

(6.27) a>y^'. 

Recall condition (I5.10p on page Uni which was derived from the closed- form solution 
of equation (16.261) and is both necessary and sufficient for (I6.26P to have a square- 
integrable solution. We conclude that, for equations with constant coefficients, (16.270 
should imply (I5.10p . but not necessarily the other way around. As a result, comparison 
of (I5.10p and (16.270 produces a lower bound on the operator norm R for the field X 
in terms of the covariance function R of the associated process X: 

^> sup ^ . 

0<t<T t 
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